oo 

o 

00 



X 



Contraction in and large time behavior for 
a system arising in chemical reactions and 
molecular motors 

Michel Chipot 0, Danielle Hilhorst 
David Kinderlehrer I, Michal Olech 

Abstract 

' We prove a contraction in property for the solutions of a nonlin- 

. ear reaction-diffusion system whose special cases include intercellular 

transport as well as reversible chemical reactions. Assuming the ex- 
istence of stationary solutions we show that the solutions stabilize as 
' t tends to infinity. Moreover, in the special case of linear reaction 

terms, we prove the existence and the uniqueness (up to a multiplica- 
' five constant) of the stationary solution. 

■ ■ Key words: weakly coupled system, molecular motor, transport, 

^ . parabolic systems, contraction property. 

; AMS subject classification: 34D23, 35K45, 35K50, 35K55, 35K57, 

^ ■ 92C37, 92C45. 

d ■ 

1 Introduction 

^ ■ We start with two specific reaction-diffusion systems. The first one 

\ describes a reversible reaction and the other one a molecular motor. 

We first consider the reversible chemical reaction (see also Bothe [1], 
■ Bothe and Hilhorst [S] , Desvillettes and Fellner [TU] and Erdi and Toth 

It involves a reaction-diffusion system of the form 

Ut = dil^u — ak{rA{u) — rsiv)^ in r2x(0, T), C , 
vt=d2/^v + (ik{rA{u)-rB{v)) in 1^ x (0,r), ^C^'^ ^"^'"^^ 
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together with the homogeneous Neumann boundary conditions, where 
di, d2, a, /?, k and T are positive constants and where is a bounded 
subset of with smooth boundary. Such systems describe, with 
a suitable choice of the functions and rs, chemical reactions for 
two mobile species. For example, functions ryi(u) = , fBiv) = v"^ 
correspond to a reversible reaction kA ^ mB. Reactions of the type 
qiAi + . . . QkAk ^ qiBi + . . . QmBm can also be described by similar 
systems with more complicated reactions terms. 

Another model problem is a system in d = 1 space dimension and 
n unknown variables ui,...,Un, n > 1, for intercellular transport, 
namely 

dui d [ dui , 



dt dx V dx 



+ ""v^^ '3r = [0, 1] X (0, T) 
a^ + Ui4 = on aQr = {0,1} X (0,r), 



where 



a-ii < 0, ttij > for all i G {1, . . . , n},i ^ j, 

" (1.2) 
aij = for all i, j S {1, . . . , n}. 

1=1 



It models transport via motor proteins in the eukaryotic cell where 
chemical energy is transduced into directed motion. A derivation of 
the system from a mass transport viewpoint is given in [7j. For an 
analysis of the steady state solutions and for further references we 
refer to [6], [12], [13], and [20]. 

In this paper we study the corresponding system in higher space 
dimension, namely 



dui - / ^ ^ I \ 

— = dw[aiVui + Ui\/^p^ ) 



+ ai(YXijrj{uj{x,t),x)] in Qt, 



(1.3a) 



where i £ {1, . . . , n}, and Ui{x, t) : Qt M+, with Qt = O x (0, T), 
ft an open bounded subset of M'^ with smooth boundary, and T some 
positive constant. We supplement this system with the Robin (no- 
flux) boundary conditions 

(Ti^+Ui^ = 0, ie{l,...,n}, on dnx{0,T), (1.3b) 
ov ou 
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where v is the outward normal vector to and the initial conditions 



ni(a;,0) = ^0,1(3;), ... , 0) = no,n(a^), x ^VL. 



(1.3c) 



We assume that the following hypotheses hold 

1. The constants Oi and Oj G M, where i € {1, . . . ,n}, are strictly 
positive; 



2. For i, j G {1, . . . , n}, Ajj < 0, \ij > if i 7^ j, X])!^;^ A^j = 0; 



3. for alH G {1, . . . , n}, the smooth functions rj are nondecreasing 
with respect to the first variable; rj(0, x) = and we assume 
that the functions V'j are smooth as well; 

4. 'Ui(.,0) = noi G C(J7), -uoi ^ 0. 

In the linear case of the molecular motors, it amounts to choosing 

rj(s, x) = s, Xij = aij and ai = 1 for all i, j G {1, . . . , n}. (1-4) 

We denote by Problem (P) the system (|1.3ap together with the 
boundary and initial conditions ()1.3bp . ()1.3cl) . and admit without 
proof that Problem (P) possesses a unique smooth and bounded so- 
lution on each time interval (0,r]. An essential idea for proving the 
existence of a solution would be to apply the Comparison principle 
Theorem 12.21 below to deduce that any solution of Problem (P) has 
to be nonnegative and bounded from above by a stationary solution. 

Finally, we note that because of the boundary conditions (|1.3bp 
the quantity 



is conserved in time. 

The organization of this paper is as follows. In Section 2 we prove 
a comparison principle for Problem (P). The main idea, which permits 
to show that Problem (P) is cooperative, is a change of functions which 
transforms the Robin boundary conditions into the homogeneous Neu- 
mann boundary conditions. In Section 3 we establish a contraction in 

property for the corresponding semigroup solution. Let us point 
out the similarity with an old result due to Crandall and Tartar [8] 
where they proved in a scalar case that in the presence of a conserva- 
tion of the integral property such as (jl.Sp , a comparison principle such 
as Theorem 12.21 is equivalent to a contraction in property such as 
the inequality (j3.4p below. As far as we know such an abstract result 
is not known in the case of systems. 

Section 4 deals with the large time behavior of the solutions. Suppos- 
ing the existence of a stationary solution, we construct a continuum of 




(1.5) 
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stationary solutions and prove that the solutions stabilize as t tends 
to infinity. Let us mention a result by Perthame [T^ who proved the 
stabilization in the case of the two component one-dimensional mo- 
lecular motor problem. Finally in Section 5, show the existence and 
uniqueness (up to a multiplicative constant) of the stationary solution 
of the molecular motor problem. 

Acknowledgment: The authors acknowledge the preliminary master 
thesis work of Aude Brisset about the corresponding two component 
system. They are grateful to the professors Piotr Biler, Stuart Hast- 
ings, Annick Lesne and Hiroshi Matano for very fruitful discussions. 

2 Comparison principle 

First, we remark that the system of equations (jl.3ap is cooperative. 
However, since nothing is known about the sign of the coefficients — 

in the Robin boundary conditions (|1.3bp . we cannot decide whether 
the Problem (P) is cooperative. This leads us to perform a change 
of variables which transforms the Robin boundary conditions into the 
homogeneous Neumann boundary conditions. 

2.1 The change of unknown functions 

Performing the change of variables 

Wiix,t)=u,{x,t)e^^^''^/''% ie{l,...,n}, (2.1) 

we deduce from (|1.3p that iv := {wi, . . . ,Wn) satisfies the parabolic 
problem 

^ = a,e'^»^^V'^«div(e-*''^'/^»Vt/;i) 
dt ^ ' 

(2.2) 

_^ ..... - , . T, 



together with the homogeneous Neumann boundary conditions 

011) ' 

-^ = 0, i£{l,...,n}, on dQ, (2.3) 
and the initial conditions 

Wi(x,0) = 'Uo,i(x)e*'^"'/"S ie{l,...,n}, x G 0. (2.4) 
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In the following, we denote by Problem P/v — the problem (j2.2p . (j2.3 
(12. 4p . To begin with we define the operators 



(2.5) 



C,iw.) = ^ - a,e'^»<^V^'div(e-^'^^V^'V^,) 

We say that {wi, . . . , u^„) is a subsolution of Problem Pjv if 

Ci{wi) ^0 in Qt, 

on dnx{0,T), (2.6) 

Wi{x,0) ^ w;j(2;,0), x G 

for alH G {1, . . . , n}. We define similarly a supersolution (ui, . . . , m„) 
of Problem P/v by the inequalities 

A(t(^j) ^0 in Qt, 

(/Til ■ 

on aJ7x(0,r), (2.7) 

Wi{x,0) ^ ?i;i(3;,0), x G fi. 

The following comparison theorem holds ([2], |21j). 

Theorem 2.1. Lei (w^i, ••• ,!«„,) and {wi, . . . ,Wn), be a sub- and 
a super - solution, respectively, for the operators Cj defined by ()2.5p 
with j G {!,..., n}, which means that (|2.6p and (|2.7p /loZc/ /or i G 
{1, . . . , n}. Then w,^ ^ uJj m Qt- Moreover, for all i G {1, . . . , n} 
stic/i i/iai w,i ^ and on {t = 0} x Q then w.^ < Wi in Qt- 



This comparison theorem immediately translates into a compar- 
ison theorem for solutions of the original Problem (P). For all i G 
{1, . . . , re}, we define the operators 

Li{ui) = {ui)t - div(fjjVuj + UiVipi) 

(2.8) 



ai(Y\ijrj{uj,x)] in Qt- 

^ .7 = 1 ^ 



The following result holds. 



Theorem 2.2. Let {u^, . . . ,u„) and {ui, . . . ,Un), be a sub - and a su- 
per-solution, respectively, for the operators Lj, defined by l\2.8\) with 
j G {1, . . . , n}. Then ^ Ui in Qt- Moreover, for all i G {1, . . . , n} 
such that Ui ^ Ui and ^Ui on {t = 0} xVt then < Ui in Qt- B 
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Next we state two immediate corollaries of Theorem 12.21 

Corollary 2.3. (uniqueness) If {u\, ... ,u\) and {u\, . . . ^u^) are solu- 
tions of Problem (P) with the same initial condition (ito,i, • • • ,UQ^n) G 
(C(II))'^, then for a// i G {1, . . . , n}, u\ = uj . M 

Corollary 2.4. (positivity) If {ui, . . . ,Un) is the solution of Problem 
(P) with the nonnegative initial condition (uo,i, • • • , uo,n) G 
then for all i G {1, . . . , n}, Ui ^ 0. Moreover, for all i G {1, . . . , n}, 
such that UQ^i ^ and UQ^i ^ 0, Ui > in Q. ■ 



3 Contraction property 

The purpose of this section is to show a contraction in (L^(r2))" pro- 
perty for the solutions of Problem (P) with the initial conditions be- 
longing to The main steps of the proof rely upon argu- 
ments due to [3] and [H]. 

We first introduce some notation. We suppose that the functions 
{u\, . . . , li^) and (uf, . . . , u^) are the solutions of Problem (P) with 
the initial conditions {u^ i, . . . , "i^o,?!) (^o,i' • • • ' '^o,n)' respectively. 
Define 

(f/l, . . . , Un) ■■= {u\ -ul,...,ui-ul). (3.1) 

Then 

iUi)t = dw{aiVUi + UiV^Ji) 

n 

+ ai'^Xij{rj{u^j{x,t),x) - rj{u'j{x,t),x)) in Qt, 

3=1 (3.2) 

cT.^ + C/.^ = on af]x(o,r), 

Ui{x,0) = Uo^i{x) for X G f], 

together with 

Uo,i = ul^i - uli, (3.3) 

for each i G {1, . . . , n}. 

Next we prove the following contraction in property. 
Theorem 3.1. For all t>0, 

— \\Ui{;t)\\mn) + ■■■ + —\\Un{-,t)\\Li(n) 

Oil Oln 

^ ;^ro,i(-)llLi(n) + • ■ ■ + ;^l|f/o,n(-)llLHf^)' (3-4) 
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where Ui and Uo^i, i G {l,...,n}, are defined by (|3.ip and (j3.3p . 

Proof Dividing each partial differential equation of (13. 2p by and 
summing them up, we obtain 



dt 



n n 

+ Aij(^rj(n](x,t),x) - rj{uj{x,t),: 

i=i j=i 

n ^ 

= V — div {aiVUi + UiVi^i) 
1=1 

i=i i=l 

= V — div (diV^/i + c/iVV'i) , 



1=1 



where we have used Hypothesis 2. 

This, together with the boundary conditions ()1.3bp . implies the 
conservation in time of the quantity 

j.±l-Ju,i.,t)i. = ,. (3.5) 

Let us look closer at the nonlinear term in (j3.2p . We can write, for 
fixed index i 



^ \ij {rj {u] {x,t),x) - rj {Uj {x,t),x)) 

n J, rj n 



i=i J=i 

Freezing the functions for i S {1, . . . , n}, A; S {1, 2}, we deduce 
that the functions C/i , . . . , C/„ satisfy a system of the form 

n 

{Ui)t = d\-v[a,VUi + t/iVV'i) + AijUj in Qt, (3.6) 

with the boundary and initial conditions 

a,^ + U,^ = Q on af7x(0,r), 

w ov (3.7) 

;7i(x,0) = C/o,i(2;), xgJ^. 
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for i G {1, . . . , n}, where Aij are functions of space and time. 
In order to make the notation more concise, we write 

Uo = {Uo,i, . . . , f^,n) ) 

U = {Ui,...,Un), 

?7± = (C/±,...,C/±), 

where = max{s,0}, s~ = max{— s,0}. By ()3.6p . ()3.7p and Corol- 
lary [2?3] we can write U in the form 

{U){x,t)=S{t)Uo{x) = {Si{t)Uo,...,Sn{t)Uo){x) 

with some operator S{t). We set 

{Wi,..., Wn) = - (C/i e^i^^V^i ,...,Un e^"'^V.^n) ^ 

and Aij = Aij e'''*'^'/'^' g-'^'i^'^'Aj _ Then, the system of equations (|3.6|) 
can be expressed in the form 

n 

(VFi)^ = cTie*'^^'/'^'div(e-'*''"'/^'VWi) + in Qt, (3.8) 

with the boundary and initial conditions 

^ = on dnx{0,T), (3.9) 
iy.(x,0) = -C/o,ie'*'*"'/"% j;GJ], (3.10) 

for i G {1, . . . , n}. 

Next we show that the solutions Wi of the problem (j3.8p - (|3.10p 
with nonpositive initial conditions are nonpositive in Q for all t G 
(0, T). To that purpose we consider the auxiliary problem 

n 

{Wi)^-Mx)dw(^Q{x)VW,)-Y,lijWj^O in Qt, (3.11) 

i=i 

on 5Ox(0,r), (3.12) 
Wi(x,0) = Wo,i(2;) < xen, (3.13) 

for i G {1, . . . ,n}. We assume that '&i{x) and 0(2^) are nonnegative 
in and that the coefficients 7ij satisfy the same assumptions as the 
coefficients Xij in Problem (P). The following result holds. 
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Lemma 3.2. Let {Wi, . . . , Wn) be a smooth and bounded solution of 
the problem ()3.1ip - (jS.lSp with nonpositive initial conditions Wo,i on 
a time interval [0,r]. Then Wi{x,t) ^ in x (0,T]. Moreover, 
for each i G {1, . . . ,n} such that Wo^i ^ and Wo^i ^ 0, Wi < in 
Qx(0,r]. 

Proof The result of Lemma 13.21 follows from the fact that the sys- 
tem (jS.lip . ()3.12p . (jS.lSp . with the inequalities {^} replaced by the 
equalities {=}, is a cooperative system. However, for the sake of com- 
pleteness, we present a proof below. We first remark that, in view of 
|21l Remark (i), p. 191], one can always satisfy the condition 

n 

jij ^ for alH G {1, . . . , n}, (3.14) 

for the matrix of coefficients (7ij)"^._j^ by performing the change of 

variables Wi = Wie~'^^ for all i G {1, . . . , n} and c > large enough. 
Thanks to the regularity of each Wi, we can apply Theorem 15, p. 191 
from [21] to conclude that - M ^ in x [0, T] for some M 
and all i G {1, . . . ,n}. In fact, we can deduce that Wi — M < in 

Q X (o,r). 

Indeed, if for some k G {l,...,n}, W^ = M in an interior point 
{x, i) gQx (0, T), then Theorem 15, p. 191 in [21] implies that Wk = 
M for all ^ t < t, which is impossible since Wk{x,0) ^ 0. If the 
maximum M of W^ is attained at a boundary point P G (90 x (0, T) 
then either there exists an open ball K C x (0, T) such that P G dK 
and Wk — M < in K, and the last part of Theorem 15, p. 191 in 
[21j contradicts the boundary inequality (I3.12p . or for all open balls 
K C ft X {Q,T) such that P G dK there exists a point (x, i) £ K such 
that Wi{x,t) = M, and we proceed as in the case before. 
Hence, there exists M > 0, such that Wi M < M inTl x [0, T] for 
all i G {1, . . . ,n}. Then we can repeat the reasoning for all M > 
until M = 0. Indeed, if this would not be the case, we find the least 
real number M > 0, with Wi ^ M ^ M in x [0,T], which leads 
again to the existence of a real number ^ M < M with the same 
property. This contradicts the fact that M was defined as the least 
such real number. ■ 

Since the functions uj, uf are bounded on x [0, T], it follows 
that the functions Wi are bounded on x [0, T] for alH G {1, . . . , n}. 
Then we are in a position to apply Lemma 13.21 with t?i(x) = ^"-l^^, 
Ciix) = (Tj e"''''/'^* and jij = Aij for i,j G {1, . . . ,n}. We deduce that 
the solutions Wi of the problem (|3.8p - (j3.1Up with nonpositive initial 
conditions are nonpositive in il. for all t G (0,T). 
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Next we remark that the above reasoning can be apphed either 
with Uq replaced by Uq or with Uq replaced by . This permits to 
show that Si{t)UQ , Si{t)UQ^ ^ and that 

Si{t)U^ > if # 0. (3.15) 
We easily compute 

n ^ 1 

— ll^^(''*)llLi(t7) - Y.—po,i{-)\\mQ) 

1=1 t=l 



1=1 * 1=1 * 

= E / -{max{5i(t)i7o+,5i(t)c7o"} 

(3.16) 

-;^min{5,(t);7o+,5,(t)c7o"}}dx-f;l / {f^+o + f^i:o} dx 
E / ;^(5.(i)f/o+ + 5.(t)C/o-)dx-f;l [ {U^o + Ur,}dx 
-2^ / -min{5i(t)?7o+,5,(t)f7o~}dx 
-2^ / — min{5i(t)?7o+,5i(t)f7o^}dx^O, 



(3.17) 

which completes the proof of ([37 



Corollary 3.3. Lei (ug . . . , Ug „), (ug i> • • • > ""o n) ^ (^(^))" '^•^ 
in Theorem Moreover, let us assume that for at least one index 
k G {1, . . . , n} t/ie difference Uqj^ — Uqi^ changes the sign. Then, the 
inequality (j3.4p is sthci for all t > 0, so that solution satisfies a strict 
contraction property. 
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4 Large time behavior of solutions 



In this section we assume the existence and uniqueness of a positive 
solution V = {vi, . . . , Vn) € {C{Vl) n C^(r2))" of the elhptic problem 

(l\\{aiVvi + ViVi}i)+aA^^\ijr,j{yj{x),x)\=^ in Q., (4.1) 

(Ji— + 't;i^— = on (90, (4.2) 

n 

V- / ^;,(x)dx = 1, (4.3) 

for i € {1, . . . , n}. 



Definition 4.1. PFe say i/iai a vector function v = {vi,...,Vn) G 

(C(r2))" is nonnegative (resp. positive) if Vi{x) ^ (resp. Vi{x) > Oj 
/or all X and all i G {1, . . . , n}. 

Next we introduce the semigroup notation for the unique solution 
of Problem (P), namely 

u{t) = T{t) no = {Ti{t) uo,...,Tn{t) uq^ , 

with the initial data uq S ((7(0))"'. The method of the proof is based 
upon an idea of Osher and Ralston |18j . It mainly exploits the contrac- 
tion properties for the nonlinear semigroup T(t) given by Theorem 
13.11 and Corollary 13.31 A similar reasoning was developed in other 
contexts by Bertsch and Hilhorst [3], Hilhorst and Hulshof [13] and 
Hilhorst and Peletier |15| . 

We suppose there exists a set C (C(r2) n (7^(0))" of positive 
stationary solutions with the following property which we denote by 
^: _ 

For each / = (/i, . . . , /n) G {C{Tl) n C2(0))" either f e Jif or there 
exists (^1, . . . , ^„) G J^, such that fi — changes the sign for at least 
one index i G {1, . . . , n}. 

One can prove that a set satisfying Property ,y exists in at 
least two cases: 

i) In the case of the system (jl.ip where the Robin boundary condi- 
tions reduce to the homogeneous Neumann boundary conditions, 



11 



the set is given by 



jr=|(a,6): a>0, b = r]^^{rA{a)) 

n 

For more details we refer to [5]- 

ii) In the case of the molecular motor with a linear n-component 
system the set is given by 

= {cv: c G M+}, 

where is a unique solution of the elliptic problem (14. ip - (14. 3p . 

Proposition 4.2. T/ie continuum zs siic/i ^/lai for each 

/ = (/i,...,/n)G (C(0)nc2(0))" 

either f £ J^f, or there exists {^i, ■ ■ ■ ,S,n) £ such that fi — 
changes the sign for at least one index i £ {1, . . . ,n}. 

Proof 

i) In the case of system (jl.ip the proof is rather obvious since the 
continuum is composed of constant pairs. 

ii) In the case of the molecular motor, let us assume that / ^ 
Then there does not exist any positive constant c such that cv = 
f. In particular, there exists an index i S {1, . . . ,n} such that 
Vi is not proportional to fi, or in other words cvi 7^ fi for all 
c > 0. Without loss of generality we can assume that the first 
coordinate has this property. Let xq € be arbitrary. Since vi 
is strictly positive in 0, we can define 

^ fijxo) 
vi{xo) ' 

so that 

(/i - CoVi){xo) = 0. 

Let 2 = |x € n : (/i — cofi)(x) = 0|. From the continu- 
ity of /i and vi, Z is closed as a subset of If there exist 
xi, X2 G Z'^, such that (/i — coVi){xi) and (/i — coUi)(x2) are 
of different signs, then the proof is complete. Now suppose that 
(/i — Co vi){x) is positive for all x G Z"^. In particular 

(/i -coVi)(x) = d > 
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for some fixed x G Z'^. Then choosing e = we see that 

2vi{x) 

{fi-{co + e)vi){x) = ^>0. 

However 

(/i - {co + e)vi){xo) < 0. 

We proceed similarly when (/i — cq f i) {x) is negative for all 
x G ■ 

In the sequel we suppose that the initial data uq = (no,i, • • • , iio,n) 
from ((7(0))" also satisfy the following property: 

There exists h € ^ such that ^ uq ^ h in i}, (4-4) 

and remark that this property is satisfied in both the cases (i) and 
(n). 

Proposition 4.3. Let uq = (uq^i, . . . , iio,n ) G (C(0))" safe/i/ i/ie 
property (j4.4p . T/ien i/ie solution (ui, . . . , u^) o/ Problem (P) is such 
that ^ ii(t) ^ /i /or all t > 0. 

Proof We remark that is a subsolution of Problem (P) and that 
/i is a supersolution, and apply Theorem 12.21 ■ 

Next we prove the main result of this section. To that purpose we 
first define the norm || • by 

" 1 

\\f 111 = I^-II^ILi{n) • 

i=i * 

Note that this norm is equivalent to the usual product norm in the 
space (L^(J7))". 

Theorem 4.4. For all nonnegative uq = (tio,i) ■ ■ ■ , iio,n) £ (C'(J^))" 
there exists f = (/i, . . . , /„) G such that 

lim II T(t)^-/ 11^=0. 

Proof 

The proof consists of several steps. To begin with we define the w-limit 
set 

uj{uo) = 1^ S i^^i^))^ '■ there exists a sequence t^ ^ oo 

as A; — > oo, such that lim ||T(tfc) — 5 |L = >, (4.5) 

k — >oo J 
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The organization of the proof is as follows. First we show that oj{uo) 
is not empty. In the second step we define the Lyapunov functional 

v(a = ||e-^||i, 

where w is a stationary solution and check that it is constant on u;{uq). 
We then deduce that u!{uo) C M', and finally prove that w(no) consists 
of exactly one function. 

Step 1. uj{uo) is not empty. 

Let e > be arbitrary. Suppose that $7' CC satisfy 

I ^ ^ 2K 



and set 

K = Y.^\\f'Mn)^ (4.6) 



n 2 



1=1 



where h has been introduced in (j4.4|) . We have already proved in 
Proposition 14.31 that T{t) uq is bounded in (L°°(r2))". Therefore there 
exist a vector function g € (L°°(J7))" and a sequence {u(tfc)} such that 

u{tk) ^ g weakly in {L\Vi)Y , (4.7) 

as tk oo. Next we deduce from [16^ Chap. Ill, Theorem 10.1] that 
there exists a positive constant C such that 

\ui{xi,t) - Ui{x2,t)\ ^ C\xi - 

for all xi, X2 G f^' and all t > 0. Therefore, it follows from the 
Ascoli-Arzela Theorem (see, e.g., [U Theorem 1.33]) that u{tk) g 
as tfc — > oo, uniformly in . We choose tQ large enough such that for 
all tfc ^ to 

H;tk)-9i-)\\i,n' (4.8) 

where || • corresponds to the norm in 0'. We deduce that, in 
view of I^Mj and dlT]) that 

||^(-,ifc)-5(-)||iMf2' ^ I' 

which together with (|4.8|) yields 

||^?(-,ifc) -5(-)||i ^ £• 



Step 2. io{uo) C 

Indeed, let g S uj{uq) and suppose ^ ^ J^. According to Proposition 



14 



14.21 we can find a steady state solution w € J^, such that at least 
one component of w — g changes the sign. Without loss of generality 
we can assume that it happens for the first component, namely that 
fi — wi changes the sign. We remark that, by the contraction property 
in Theorem 13. 11 the functional 

v(a = ik"'-4i 

is a Lyapunov functional for Problem (P), where ^ G Next 
we describe some of its properties. 

Property (a) The functional V is constant onu>{u()). 

Since T{t)w = w and T{t) has the contraction property (j3.4p . the 

functional V is nonincreasing in time along the trajectory T{t)uQ, 
which yields 

V{T{t)uo) = \\T{t)uo-w\\^ 

= \\T{t) uq — T{t) w\\^ =^ ll^^o — ^lli < oo . 

Thus there exists a finite limit V* of V(T(t)uQj as t oo. Let 
hi, /i2 € ^{uq). We can find a sequence ^ oo as /c — > oo, such that 

\\'^{t2k)uo -hi\\^ ^ and ||T(i2fc+i) - ^2||i ^ 0, 

as k tends to oo. It follows that V(/ii) = V(/i2) = V*. 

Property (b) The uj-limit set uj{uo) is invariant with respect to the 
semigroup T{t) , namely if /i € uj{uq), then for all t > also T{t) h G 
io{uo). 

Let the sequence — > cxd as A; — > oo be such that uq — h\\i 0. 

Prom the contraction property (|3.4|) 

\\T{tk + t)uo-T{t) h\\^ = \\T{t)T{tk)uo-T{t)h\\^ 
^ \\T{tk)uo-h\\^. 

Since the last term above tends to as /c tends to oo this shows that 
T{t)h£ Lu{uo). 

Now, remember that g G uj{uo) is such that g ^ and w G Jif is 
such that the first component of w — g changes the sign in 0,. Then, 
Corollary 13.31 yields 

V{T{t)g) = \\T{t)g-w\\^ 

= \\Tit)g-T{t)w\\^ < \\g-w\\i = Vig), 
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for all t > 0, which contradicts Property (a). Therefore g € Jif. 

Step 3. The set uj{uq) contains only one element. 
Suppose that 51,92 G "^(^^o)- Then we can find two sequences tk, Sk 
tending to 00 as —s- 00, such that s^ ^ tk and ito— ||'^('Sfc) uq — 

^ as tk 00. Since uj{uo) C J^, it follows that 

Wgi - 52II1 ^ ll'^(ifc) uo - gi\\^ + \\T{tk) no - 52||i 
= \\T{tk)uo - gi\\i + - Sk)T{sk)uo - T{tk - Sfc)52||i 

^ ||T(tfc) no - + ^T{sk) no - 52II1 , 

which tends to as /c ^ 00. ■ 



5 Stationary solutions for the linear 
molecular motor problem 

In this section we show the existence and the uniqueness (up to a mul- 
tiplicative constant) of the classical stationary solution of the problem 
for the molecular motor. We suppose that is an open bounded sub- 
set of M'^ with smooth boundary 50. 
We consider the linear system 

n 

div{aiVvi{x) + Vi{x)V'tl)i{x)) +'^^aijVj{x) = in Vt, (5.1) 

i=i 

where i G {1, . . . ,n}, n > 1. The system (|5.ip is supplemented with 
the Robin boundary conditions 

a.$^+n,5^ = on d^, (5.2) 

where i G {1, . . . , n}. Thus, the problem can be written as 

Av = 0, 

with a linear operator ^ in a suitable Banach space X of functions 
on 0, to be made precise later. Moreover, we impose the integral 
constraint 

n „ 

/ Vi{x)dx = 1. (5.3) 
The adjoint problem A*(p = to (15. ip . in a dual space X*, is now 



aiAifi - V^i • V(pi + ^ CLjiipj = 0, in 0, (5.4) 
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with the Neumann boundary conditions for each i = 1, . . . , n 

^ = on on. (5.5) 
Since X]j=i '^ji — 0) the problem ()5.4p has the obvious solution 

if={ipi,...,ipn) = {l,...,l)- (5.6) 

We are going to apply the Krein-Rutman theorem on the first eigen- 
values and eigenvectors of positive operators, and this will permit us 
to conclude that the problem ()5.1l) - ()5.2p has a one-dimensional space 
of solutions. Therefore, under the additional constraint (|5.3p . the orig- 
inal problem ()5.ip - ()5.2p has a unique solution. 

Perthame and Souganidis sketched this argument for n > 1 and d = 1 
in [20]. 

Theorem 5.1. Under the assumption Yl^=i ^ji ~ 0, there exists a 
unique smooth solution v of the system ll5.1\) - [^73\) . 

Before proving Theorem 15. II we recall some basic definitions as well 
as the Krein-Rutman theorem from [9l Ch. VIII, p. 188-191]. 

Definition 5.2 (Reproducing cone). We say that a closed set K in 
X is a cone, if it possesses the following properties: 

i) oeK, 

a) u, V ^ K =^ au + f3v ^ K , for all a, (3^0, 
Hi) V K and —v E K =^ v = 0. 
A cone K (Z X is said to be reproducing if X = K — K = {^ki — k2 '■ 

h, K}. 

Definition 5.3 (Dual cone). If K is a cone in X, then the set K* C 
X* is said to he a dual cone if 

{f*,v)>0, 

for every v (z K. 

Definition 5.4 (Strict positivity). Let B be a linear operator on X. 
Then B is said to be strongly positive if Bv € K° for all v G K such 
that f / 0. 

Theorem 5.5. Let K be a reproducing cone in a Banach space X , 
with nonempty interior K° ^ 0, and let B be a strongly positive com- 
pact operator on K in a sense of Definition \5.4\ Then the spectral 
radius of B, r{B), is a simple eigenvalue of B and B* , and their as- 
sociated eigenvectors belong to K° and {K*)° . More precisely, there 
exists a unique associated eigenvector in K° (resp. (K*)") of norm 1. 
Furthermore, all other eigenvalues are strictly less in absolute value 
than r{B). 
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Proof We will apply Theorem 15.51 to the space Af = ((7(0))"' C 
(L^(r2)) endowed with the usual supremum norm, and the operators 

B= {XI- Ay^ : X ^ X, 
B* = {\I-A*Y^ : X* X\ 

where A > is a strictly positive real number to be fixed later. 
Let 

K = ^ X : Ui{x) ^ for each x G 0, i = 1, . . . , n}. 
We remark that is a reproducing cone, with nonempty interior 
K° = [u^X: in£'Uj(x) > 0, -j = 1, . . . 

From the standard theory |17l Theorem 2.1 and Theorem 3.1, Ch. 
7] for elliptic partial differential linear systems, the boundary value 
problem 

n 

Oil^Lpi - Vlpi ■ Vifi + ^ ajiifj - Xipi = fi in Q, (5.7) 
1=1 

with the homogeneous Neumann conditions (|5.5|) on dfl, for A = A > 
sufficiently large, has a solution = (ipi, . . . ,ipn) € X for each 
/ = (/i, . . . , fn) G X. Moreover, if fi{x) ^ for each i = 1, . . . ,n, 
and X G fi, then ifi{x) ^ (in fact, (pi{x) > in fi), which is a 
consequence of the maximum principle (cf. also Example 3 on p. 196- 
197 in [9]). Thus, the operator B* = [XI — A*) is a strongly positive 
and compact operator, and by Theorem 15. 5^ the largest eigenvalue // 
of B and B* is simple. 
Since 

n 

-aiAipi + V^jJi ■ Vifi - ^ ajiipj + Xipi = Xipi in Q, 

^'^i n fin 
— — =0 on oil, 

ou 

for all i G {1, . . . , n}, with if = (991, . . . , ipn) = (1, . . . , 1), and since 

(1, . . . , 1) G {K*Y, it follows that i = r({XI - A*)~^\ is a simple 

~ 1 

eigenvalue of the operator [XI — A*) . Applying again Theorem 15.51 

we deduce that — is the largest eigenvalue of the operator (A/ — A) ^ 
X 

and that it is simple, and that there exists v G K" C X such that 
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which is equivalent to 

Av = 0. 

This proves the existence of the solution of the problem (j5.ip - (|5.3 



References 

[1] R. A. Adams, Sobolev Spaces, Pure and Applied Mathematics, 
vol. 65, Academic Press, New York-London, 1975. 

[2] N. D. Alikakos, P. Hess, H. Matano, Discrete order preserving 
semigroups and stability for periodic parabolic differential equa- 
tions, J. Differential Equations, 1989, 2, 82, 322-341. 

[3] M. Bertsch and D. Hilhorst, A density dependent diffusion 
equation in population dynamics: stabilization to equilibrium, 
SIAM J. Math. Anal, 1986, 4, 17, 863-883. 

[4] D. Bothe, Instantaneous limits of reversible chemical reactions 
in presence of macroscopic convection, J. Differential Equations, 

2003, 1, 193, 27-48. 

[5] D. Bothe and D. Hilhorst, A reaction-diffusion system with fast 
reversible reaction, J. Math. Anal. AppL, 2003, 1, 268, 125-135. 

[6] M. Chipot, S. Hastings, and D. Kinderlehrer, Transport in a 
molecular motor system, M2AN Math. Model. Numer. Anal., 

2004, 38, 6, 1011-1034. 

[7] M. Chipot, D. Kinderlehrer, and M. Kowalczyk, 2004 A varia- 
tional principle for molecular motors, Meccanica, 38, 505-518 

[8] M. G. Crandall and L. Tartar, Some relations between nonex- 
pansive and order preserving mappings, Proc. Amer. Math. Soc, 
1980, 3, 78, 385-390. 

[9] R. Dautray, J.-L. Lions, Mathematical Analysis and Numerical 
Methods for Science and Technology, Volume 3: Spectral Theory 
and Applications, Springer- Verlag, Berlin, 1990. 

[10] L. Desvillettes and K. Fellner, Exponential decay toward equi- 
librium via entropy methods for reaction-diffusion equations, 
J. Math. Anal. AppL, 2006, 1, 319, 157-176. 

[11] P. Erdi and J. Toth, Mathematical Models of Chemical Reac- 
tions, Nonlinear Science: Theory and Applications, Princeton 
University Press, Princeton, NJ, 1989. 



19 



[12] Hastings, S., Kinder lehrer, D. and McLeod, J.B. Diffusion me- 
diated transport in multiple state systems, An elliptic-parabolic 
problem in combustion theory: convergence to travelling waves, 
SI AM J. Math. Anal, 2007, 39.4, 1208-1230. 

[13] Hastings, S., Kinderlehrer, D. and McLeod, J.B. Diffusion medi- 
ated transport with a look at motor proteins, Recent advances in 
nonlinear analysis, (Chipot, M., Lin, C-S,., and Tsai, D-H., eds) 
(2008) World Scientific, 95 - 112. 

[14] D. Hilhorst and J. Hulshof, An elliptic-parabolic problem in com- 
bustion theory: convergence to travelling waves, Nonlinear Anal- 
ysis TMA, 1991, 6, 17, 519-546. 

[15] D. Hilhorst and M. A. Peletier, Convergence to travelling waves 
in a reaction-diffusion system arising in contaminant transport, 
J. Differential Equations, 2000, 1, 163, 89-112. 

[16] O. A. Ladyzenskaja and V. A. Solonnikov and N. N. Ural'ceva, 
Linear and Quasilinear Equations of Parabolic Type, Translations 
of Mathematical Monographs, Vol. 23, American Mathematical 
Society, Providence, R.I., 1967. 

[17] O. A. Ladyzenskaja and N. N. Ural'ceva, Linear and Quasilinear 
Elliptic Equations, Academic Press, New York, 1968. 

[18] S. Osher and J. Ralston, stability of travelling waves with 
applications to convective porous media flow, Comm. Pure Appl. 
Math., 1982, 6, 35, 737-749. 

[19] B. Perthame, The general relative entropy principle — 
applications in Perron-Frobenius and Floquet theories and 
a parabolic system for biomotors. Rend. Accad. Naz. Sci. XL 
Mem. Mat. Appl. (5), 2005, 29, 307-325. 

[20] B. Perthame, P. E. Souganidis, Asymmetric potentials and motor 
effect: a large deviation approach, 1-16, Arch. Rat. Mech. Anal., 
to appear, 2008. 

[21] M. H. Protter and H. F. Weinberger, Maximum Principles in 
Differential Equations, Springer- Verlag, New York, 1984. 



20 



